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ABSTRACT 


Two unsteady flows donunated by the occurrence of separation are simulated 
through the use of the discrete vortex model. The first of these is a sinusoidally 
oscillating flow about a circular cylinder at a Keulegan-Carpenter number of K= 10. 
The vortex model has been combined with the boundary layer calculations and the 
positions of the separation and stagnation points, the evolution of the wake, the 
velocity and pressure distributions, and the instantaneous forces have been calculated 
and compared, whenever possible. with those obtained experimentally. The model has 
successfully siniulated the occurrence of the transverse half Karman Vortex Streersmaad 
calculated positions of the vortices were found to be in good agreement with those 
obtained experimentally. The measured and calculated in-line forces and the 
differential pressure distributions showed reasonably good agreement. 

The second simulation dealt with a rapidly decelerating flow about a two- 
dimensional sharp-edged camber. An extensive study of the velocity field in the 
vicinity of the singular points led to the development of a novel method for the 
introduction of vorticity at variable time intervals. The measured and calculated 
characteristics of the flow, such as the evolution of the wake and the forces acting on 
the camber, were found to be in excellent agreement. Furthermore, the Sinimlagiog 
provided a plausible explanation for the cause of parachute collapse. a phenomenon 


which gave impetus to the numerical and physical experiments described herein. 
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I. DISCRETE VORTEX ANALYSIS OF SINUSOIDALLY OSCILLATING 
FLOW 
ABOUT CIRCULAR CYLINDERS 


A. INTRODUCTION 

The separated steady and unsteady flows about bluff bodies have been almost 
completely unvielding to both analysis and numerical simulation for a number of 
mathematical reasons and fundamental fluid dynamic phenomena. Separation gives 
rise to the formation of free shear layers which roll up into vortex rings or counter- 
rotating vortices. They, in turn, interact with each other, with the counter-sign vorticity 
generated at the base of the body, and with the motion of often unknown separation 
points. The wake becomes unsteady even for a steady ambient flow and the problem 
of the determination of the characteristics of the wake becomes coupled to the 
conditions prevailing upstream of the separation points. Evidently. viscositv modifies 
radically the inviscid flow, which, in this case, cannot serve even aS gue 
approximation to the actual flow. The boundary layer equations are not applicable 
bevond the separation points and are, therefore, of limited use in bluff-body flow 
problems. 

Fage and Johansen’s pioneering experimental work (1928). Gerrard's (1966) 
Vortex formation model, and Roshko’s (1954) numerous contributions, followed by a 
large number of important papers, have provided extremely useful insights into the 
mechanism of vortex shedding. It became clear that a two-dimensional bodv immersed 
in a two-dimensional steady tlow does not give rise to a two-dimensional steadv wake 
and oniy a fraction {about 60% for a circular cylinder) of the original-cireuianeme 
survives the vortex formation. It also became clear that bluff-body flows exhibiting 
separation, turbulence, and time-dependence are almost completely unvielding to both 
analysis and simulation even if the ambient flow is assumed to be time invariant. 

Many flows of practical interest are unsteady, i.e., the characteristics Omenae 
ambient flow are time-dependent. In the past twenty vears or so a large number of 
theoretical and experimental studies have been carried out. These dealt primarily with 
unseparated laminar flows, the early stages of impulsively started flow over plates and 
cvlinders (numerical and experimental studies), and oscillating flows with zero or non- 


zero mean flow (on an infinite flat plate and over a cylinder with streaming flow, all 


i4 


under laminar flow conditions) for the purpose of studving the effects of flow 
unsteadiness on the transition mechanisn. and turbulence development (see e.g. 
Bradbury et al. 1982). Very little has been attempted either theoretically or 
experimentally to analyze the wake-boundary-layer interaction in time-dependent flows 
(.e., with unsteady ambient flow). 

he subject of separated time-dependent flow at large Reynolds numbers 1s lesser 
developed but of greater practical importance (particularly to marine related topics) 
relative to other classical coniponent disciplines of fluid mechanics. 

A number of unsteady flow machines and their use in the investigation of 
unsteady turbulent boundary layers have been reviewed and documented by Carr 
(1981). These included flat plate, diffuser, pipe. airfoil, and cascade flows. The results 
have shown that (1) the time-averaged miean velocity profile 1s almost always the same 
as the velocity profile that would occur in a steady flow having an equivalent mean 
external flow velocity; (11) the turbulent structure in the oscillating flow 1s not changed 
from the equivalent steady-state counter part; and (in) the unsteady effects are often 
confined to a thin layer near the wall, while the outer region of the boundary laver 1s 
not strongly affected. These conclusions, apparently valid for unsteady turbulent 
boundary layer flows, are not applicable to unsteady separated, turbulent, bluff-bodv 
flows. 

The separated unsteady flow situations involving wake return, as in the case of a 
sinusoidally oscillating flow about a cylinder, or wake retardation, as in the case of a 
decelerating parachute, are an order of magnitude more complex. 

In steadv flow the position of the separation points is nearly stationary, except 
for small excursions of about 3 degrees (on a circular cylinder). Furthermore, the 
Mmerierence between the vortices and the body is confined mostly to the vortex 
formation region. 

omosciianimestiows the met eflect of the shed vortices is twofold. Firstly, their 
return to the body dramatically affects the boundary laver, outer flow, pressure 
distribution, and the generation and survival rate of the new vorticity. Secondly, they 
not only give rise to additional separation points (during the early stages of the flow 
reversal) but also strongly affect the motion of the primary separation points. These 
effects are further compounded by the diffusion and decay of vortices and by the three- 
dimensional nature of the flow (all of which give rise to cycle-to-cycle variations, 


Mmemicrous flow miodes, etc.). The stronger and better correlated the returning vortices, 


the sharper and more pronounced the changes are in the pressure distribution on the 
body and in the integrated quantities such as lift, drag, and inertia coefficients. 
Nevertheless, the increased correlation does not entirely eliminate the consequencesson 
the stochastic variations in the motion of vortices. 

in periodic flow, the mobile separation points (when they are not fixed by sharp 
edges}, undergo large excursions (as much as 120 degrees during a siven @veleuan 
oscillating flow over a circular cylinder). This experimental fact renders the treatment 
of boundarv lavers on bluff-bodies subjected to periodic wake return extremely difficult, 
particularly when the state of the boundary laver changes during a given cycle. 
Furthermore, the classical criterion of separation for steady flow, i.e., the vanishing of 
skin friction on the body, is no longer valid for unsteady flow. According to the MRS 
criterion (Moore 1958, Rott 1956, and Sears 1972), it is the simultaneous Vanishimemem 
the shear and velocity at a point within the boundary layer that determines the 
separation point. Furthermore, the time rate of change of circulation is no longer 
een dP dt=0.5U.- as in steady flow, where L. is the outer flow velocity at 
separation, but by (0.5 U- = U8.) where 08). is the speed of the separation point. It is 
clear from the foregoing that there is little hope of devising a satisfactorv theoretical 
model before something is understood of the unsteady processes associated with the 
formation and reversal of the wake, spanwise coherence, the sensitive dependence of 
the motion of vortices on small changes in the previous conditions and on the nature 


of transition in oscillating flow about smooth and rough cvlinders. 


B. NUMERICAL SOLUTION OF THE NAVIER-STOKES EQUATIONS 

For steadv ambient flow about bluff bodies (mostly cylindrical and airfoils), the 
numerical studies based on the use of the steady or unsteady form of the Navier-Stokes 
equations and some suitable spatial and temporal differencing schemes are linuted, out 
of necessity, to relatively low Reynolds numbers (less than about 1,000) (see eiey 
Lecointe & Piquet (1984) for a finite difference solution and Gresho et al. (1984) for a 
finite element solution). The major obstacles to the application of either thesia 
difference or the finite element methods to higher Reynolds number laminar flow are 
stabilitv, computation ume, treatment of the boundary conditions, and accuracy. Even 
though many differencing schemes have been developed to overcome the instability 
problem (Roache 1976), maintaining stability continues to be a problem with 
increasing Reynolds number. The truncation errors decrease the apparent Reynolds 


number by introducing an unknown artificial viscosity. Even if the problems 
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associated with stability and truncation errors were to be resolved, the attempt to 
obtain solutions of the Navier-Stokes equations at higher Reynolds numbers are 
limited by a fundamental fluid dvnamic phenomenon: the stability of the flow itself. 
When the flow becomes turbulent either in the wake and/or in the boundary layers, 
one needs a closure model for turbulence to solve the Reynolds equations for a time- 
dependent, three-dimensional, separated, turbulent flow (even if the ambient flow is 
smooth and the bluff body is two dimensional). Clearly, the roots of the most serious 
problem in the solution of the Navier-Stokes equations are buried in the physics of 
turbulence. The stability and truncation-error problems associated with the 
differencing schemes may be resolved but the problem of turbulence appears to 
transcend all efforts. 

The numerical solution of unsteady incompressible Navier-Stokes equations in 
their vorticity-stream-function formulation has been investigated by numerous 
feseatcners through the use of various finite-difference techniques. These studies 
concern miostly the separated flow about circular cylinder and prisms at relatively low 
Feynolds numibers (see e.g., Davis & Vioore 1982). 

It appears that the existing numerical methods cannot yet treat the high 
Reynolds number flows with sufficient accuracy for a number of reasons. The finite 
difference schemes require a very fine grid, a turbulent model, and a very large 
computer memory. It seems that the modelling of the turbulent stresses in the wake, 
particularly in time-dependent flows, will be the major source of difficulty in all future 
calculations. Whether or not it will ever be practical to apply the finite difference and 
finite element methods to high Reynolds number flows is unknown. The inherent 
difficulties are certainly significant enough to warrant exploring other solution 


methods. 


C. VORTEX METHODS 
1. Discrete Vortex Analysis 

Certain separated time-dependent flows may be simulated through the use of 
Meemciscrete vortex model (DV Vi) (see e.g. Chorin 1973; Sarpkava 1975). The free 
shear lavers which emanate from the sides of the body are represented by an assembly 
of discrete vortices. It has not yet been proven that a continuous vortex sheet may be 
so discretized. Thus, attention is given here to large scale vortex structures rather than 
to small-scale instabilities resulting from the vortex interactions. Furthermore, only the 
most important features of the method, as it is applied to sinusoidally oscillating flow 


about a smooth cylinder are described. 
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The previous attempts to apply the DVM to oscillating flow had either gross 
simplifications or met with various difficulties. Ward and Dalton (1969) considered 
only the symmetric flow situation with fixed separation points. Stansby (1977, 1979, 
1981) fixed the separation points at + 90 degrees and used the velocity of the nascent 
vortex rather than the velocity at the separation point to calculate the vortex strengths. 
This resulted in significantly less vorticity input and prevented the returning vortices 
from interacting freely with the boundary layers and separation points. Subsequently, 
Stansby and Dixon (1983) used a Lagrangian vortex scheme and replaced the body 
surface by a polygon of line segments. The strengths of the segments of vortex sheet 
needed to establish the zero-velocity condition on the surface were determined frorn the 
inverse of an influence matrix and the tangential velocities just inside the summed 
Subsequently, the segments were replaced by one or more point vortices. Stansby's 
calculations for Keulegan-Carpenter number K = 10, (K = U_T/D), failed to predict 
the transverse Vortex Street obServed experimenting 

Sawaragi and Nakamura (1979) determined the separation points using 
Schlichting’s (1932) periodic boundary laver theory (valid only for K< <1). Thev have 
not used the Kutta condition and incorrectly included an image vortex at the center of 
the cvlinder. Finally, the calculations were performed for only three-quarters of a 
cycle, hardly enough time for the transient flow to develop into quasi-steady state. 
Kudo (1979, 1981) investigated the sinusoidally oscillating flow about a flat plate 
normal to the flow. The wake was assumed to remain symmetrical. Kudo’s model 
used a Kutta condition, combined with a highly complicated force-and momentuim-free 
nascent vortex placement scheme. Ikeda and Himeno (1981) studied the oscillating 
flow about a cvlinder and a Lewis form. Separation points were assumed to be given 
by Schlichting’s (1932) solution. As with Sawaragi and Nakamura (1979), they have 
incorrectly retained the image vortices at the center of the cylinder. 

The discrete vortex model, as used in the present investigation, is relatively 
free from the arbitrary assumptions and inconsistencies noted above. Ilowever, the 
inclusion of the effect of turbulence in the boundary lavers and in the wake remains 
unresolved. As it will be noted shortly. the determination of the separation points 
requires the use of a separation criteria based on a laminar or turbulent flow 


separauon. 
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2. Formulation of the Problem 


The complex velocity potential may be written as 


ee 
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where U(t) = sin (2mt’T), F, and z, are the strength and position of the n-th vortex. 
The velocities are normalized by U_ and the distances by c. The complex velocity 1s 
even bred idz = —u + Iv. 

The instantaneous force acting on the cylinder may be calculated either 
through the integration of pressure or through the use of the rate of change of impulse. 


It is relatively easy to show that the pressure on the cylinder is given by 
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In terms of dimensional quantities, the use of the rate of change of impulse 
yields, 
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where D and L represent, respectively, the in-line and transverse forces. 


3. Specific Details of Creation and Convection of the Vortices 
The solution procedure employed was as follows: 

(1) The positions of the stagnations points at the upstream and downstream faces 
of the cylinder are calculated. For this purpose, the points at which the 
velocity 1s zero and changes its sign are located, starting from the most recent 
Stagnation points. When the ambient flow reverses its direction (at the 
deginning of each cycle) the stagnation points switch their positions (the 


upstream one becomes downstream and Vice versa). 


I? 


(2) The positions of the primary and secondary separation points are calculated 


through the use of one of the following methods: 


(a) 


(b) 


Pohlhausen’s Method: The velocity is calculated at one degree intervals 
along the cvlinder, starting at the stagnation points. Then the positions 
of the separation points are determined through the use of Pohlhausen’s 
technique (for details see e.g., Schlichting 1932). When the Womeiees 
returning to the cylinder cause irregularities in the velocity distribution 
along the cvlinder (because of the disproportionatelv large influence of 
some of the point vortices), the Pohlhausen’s method fails to predict a 
separation point. Under these circumstances, the maximum velocity 
Criletianals vised: 

Maximum Velocity Criterion: The velocity distribution is calculated 
along the cylinder and the position of the maximum absolute velocity on 
each side of the cylinder 1s determined. Then the separation points are 
located at points where the tangential velocity is a certain time-averaged 
fraction of the maximum velocity. Further details of this fraction will be 
discussed later. 

Absolute Maximum Pressure Criterion: When the difference between 
the primary separation angles calculated at umes t, and t__,. either 
through the use of Pohlhausen’s technique or through the use of the 
maximum Velocity criterion. is larger than 8 degrees, then the positions 
of the maximum absolute pressure are used to calculate the separation 
points. For this purpose, all velocitv peaks are calculated through the 
use of the tangential velocity distribution. Then the maximum velocity 
at which the maximum absolute pressure occurs 1s determined. [hen the 
positions of the separation points are calculated using the maximum 
velocity criterion. In other words, the pressure distribution 18 Useqmem 
locate the position of the maximum velocity nearest the true separation 


point. 


The foregoing, relatively time consuming, steps could have been elinunated 


through the use of a suitable numerical filter so as to remove the secondary 


oscillations imposed on the velocity distribution by a few point vortices in the 


vicinity of the cylinder. It was realized that the use of such a filter willemor 


only consume more computer time but will also introduce uncertainties into 


the calculations. Thus, its use was disregarded. 
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(4) 


(0) 


The velocity U, is calculated at each separation point and the ratio LU. is 
determined. Then the cumulative average of this ratio is calculated for use in 
connection with the maximum Velocity criterion to determine the position of 
the separation points whenever the Pohlhausen criterion fails. 


The strength of the primary nascent vortex is calculated from, 


op = (0.5 1U,| U,—U,08,)-At (1.4) 


where U.c8. accounts for the motion of the separation point relative to the 
flow velocity prevailing at the separation point. The strength of the secondary 


vortex 1s calculated from. 


[T..= 0.5 |Uj U At (1.5) 


In other words, the effect of the relative motion of the secondary separation 


point on I’, is ignored. 


The nascent vortices are placed at a distance € along the radial lines passing 
through the separation points, ie., atz, = (1+) exp(i @.). The value of ¢ 
which satisfies the Kutta condition is given by (see e.g., Sarpkava & Shoaff 
1979) 
reat 2 UL) 
Ce — Tae (1.6) 
LS ( He j 2 mL.) 
The velocity at any point r<r* 1s calculated using the exact solution of the 


Navier-Stokes equations for a single rectilinear viscous vortex (Lamb 1932). 
1e., -utiv=(P;2nr) [1 — exp( — 7 dvt)] e!? and at any point r>r* using the 
point vortex relation {(-ut+iv=dw//dz). The radius r is measured from the 
center of the vortex and 0 is taken plus in the CCW direction from the x axis. 
The core radius r* at which the tangential velocitv is maximum is given by 
r*=./1.26(vt) where t is the age of the vortex since its inception. Then the 
vortices are convected through the use of a simple Eulerian scheme for a time 
interval At = 0.125¢;U' T 


(7) The vortices which come closer than a distance of 0.05c to the cylinder surface 
as Well as those which have a strength smaller than 0.005 are removed from 
the field. 

(S) Oppositely signed vortices are combined when their separation distance is less 
than 0.1. Furthermore, the secondary vortices are combined at suitable 
intervals with the primary vortices to reduce the number of the vortices and 
hence the computation time. 

(9) The pressure and tangential velocity distributions and the force acting on the 
cvlinder (both through the integration of pressure and the rate of change of 
impulse, given by Eq: (1.3)) are calculated 

(10) At the end of the foregoing calculations, plots of the vortex positions, 
tangential velocity, total pressure, and stagnation and separation points are 
produced. 

Two additional details of the model need to be described: The introduction of 
asymimetry and the circulation dissipation. Anv flow started impulsively from rest 
remains svnimetrical (at least in its numerical simulations) if not disturbed for a short 
time period at or near the beginning of its inception. Sarpkaya and Shoaff (1979) have 
investigated the methods of asymmetry introduction and found that the displacement 
of the vortices in one shear layer for a short time interval 1s much more suitable than 
any of the methods previously used. In the present calculations, a similar method was 
used. The only difference was that asymmietry was introduced during roughly the first 
quarter cycle of flow (0.1 <t:T <0.3) in order to allow the asymmetry to take effecutiima 
reasonable time period. 

The previous investigations by Sarpkaya & Shoaff (1979) and others (e.g... 
Kiva & Arie 1977) have shown that the incorporation of dissipation into a discrete 
vortex model resulted in a reduction of lift and drag force magnitudes, but did not 
significantly affect the flow kinematics. In the present work, sample calculations with 
or without dissipation have shown also that the results differ only in the magnitude of 
the force coefficients and that the kinematics of the flow is not affected by dissipation. 
It is in view of this fact that in the results presented herein it was decided to avoid a 
relatively subjective dissipation mechanism. 

4. Discussion of Results 
The numerical calculations were carried out for a Keulegan-Carpenter number 


of K = 10. This value of K was chosen primarily because of the fact that some of the 


I 
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niost important phenomena take place in the range 8< K< [5 (Sarpkaya 1985: Bearman 
1985: Williamson 1985). The most important of these phenomena is the occurrence of 
a transverse half Karman vortex street on one side of the cylinder. 

A series of flow visualization experiments at K = 10 was conducted in a water 
Paplemericures Ileana I2 show a time sequence from a representative run. In these 
experiments the cylinder (D = 1.5 inches) was oscillated in the water table using a 
period of T = 3 seconds (Re = 12.800 and Re/K = B = 1280). The end of the 
plexiglass cylinder was within 1,16 inch of the bottom of the water table, so that the 
end effects were minimal. It can be seen that the vortex shedding 1s indeed on the 
same side (1.e., left) of the cvlinder. The numerical simulation of this phenomenon 
ipeomeme tie use of the discrete vortex model became a challenge to numerous workers 
(Sawaragi & Nakamura 1979; Ikeda & Himeno 1981; Stansby 1977, 1979, 1981; and 
Stansbv & Dixon 1982, 1983). However, none of these investigators has been able to 
simulate the half Karman vortex street for a number of reasons described in the 
Introduction. The present simulation has removed the shortcomings of the previous 
analyses, discovered physically and theoretically realistic methods to deal with the 
Separation points. eliminated the ad-hoc assumptions, and produced the results 
presented herein. 

Figures 1.3 through 1.7 show, at times T* = t/T = 0.2, 0.4875, 0.6063, 
ieye25, and 0.8625 the position of the discrete vortices, the tangential velocity 
distribution and the total pressure distribution. The arrow in each circle shows the 
direction and the magnitude of the ambient velocity. The vortices of opposite 
circulation are shown with different symbols. These figures will be used to familiarize 
the reader with the essential elements of the analysis. Subsequently, the evolution of 
the transverse vortex street, characteristics of the shear layers, and the normalized drag 
and lift forces will be presented. 

As the ambient velocity starts from rest and nears its maximum value (Fig. 
1.3), the vortices grow nearly symmetrically on the downstream side of the cylinder. 
reminiscent of an impulsively started flow about a cylinder. This is a consequence of 
the spiralling of the primary shear layers (SL-1 and SL-2), emanating from the primary 
Separation points (SP-1 and SP-2). The primary vortices give rise to secondary 
boundary layers on the downstream side of the cvlinder. These, in turn, separate and 
give rise to two additional shear lavers (called the secondary shear lavers, denoted by 


Sueeseand SL-4). As time increases and the ambient velocity decreases (Fig. 1.4), the 
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first vortex sheds from the upper right hand side of the cylinder. This leads to the 
rapid growth of the vortex in the lower right hand side. The tangential velocity around 
the cylinder is now primarily due to that induced by the vortices. The spikes in the 
velocity and pressure distributions are a consequence of the disproportionately large 
influence of a few vortices in the vicinity of the cylinder. Thev are most prevalemmenem 
times of very small ambient velocity. 

For T* > 0.5, the ambient flow reverses its direction, the flow Separaresmem 
the left side of the cylinder, and the previously shed vortices are convected towardsslem 
(Figs. 1.5 and 1.6) partly due to the ambient flow and partly due to the niutual 
induction of all the vortices and their images. The reversal of the flow has two major 
consequences. First, the convection of the previously shed vortex over the shoulemsen 
the cvlinder precipitates earlier separation and establishes, by its sense of rotation, a 
preferred position for the next donunant vortex. Second, the reversal of flow gives rise 
to additional primary and secondary separation points (Fig. 1.6). However, at a given 
time there are at most four shear layers, i.e., some of the separation points disappear 
and some new ones come into existence as the flow reverses. Subsequcntly ee 
vortical structure at the upper left side of the cvlinder moves further downstream and 
the vortex on the lower left side of the cylinder (Fig. 1.7) grows rapidly (in a manner 
simular to that shown in Fig. 1.4). It begins to move towards the top of the cylinder, 
partly due to the mutual induction velocity of the aforementioned vortical structure 
and partly due to that of its image. The events just described miore or less repeat 
thiemselves in the subsequent cycles (see Figs. A.1 through A.17 in Appendix Aj syiGne 
or less, because no two cycles can be expected to be exactly alike Gue Tomtiemeus 
increasing numiber of vortices or, in other words, due to the transient nature of the 
flow (for an experimental confirmation of this fact see Sarpkaya 1986). Many more 
cycles of calculations will have to be carried out in order to reach a quasi-periodic 
state. However, this is nearly impossible and hardly necessary, partly because of the 
computer-time linutations and partly because the calculations over a period of three 
cvcles are more than sufficient to delineate the fundamental characteristics Gf thesiigms 
The reason for the latter is provided by the mechanism of vortex shedding during each 
half cvcle of flow. As noted in the course of the discussion of Figs. 1.3 through 1.7, 
the flow in each half cycle behaves as if it started more Or 16ss) anew tics ima 
convected vortex over the shoulder of the cvlinder precipitating earlier separation and 


dictating the side from which the next dominant vortex is to be shed. In numerical 
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Figure 1.3 Position of vortices, velocity and pressure distributions at T*=0.200. 
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Figure 1.4 Position of vortices, velocity and pressure distributions at T* = 0.4875. 
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Figure 1.5 Position of vortices, velocity and pressure distributions at T* = 0.6063. 
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Figure 1.6 


Position of vortices, velocity and pressure distributions at T*=0,7625. 
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Figure 1.7 Position of vortices, velocity and pressure distributions at T*= 0.8625. 
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calculations the vortices do not dissipate and, therefore, begin to cover the entire flow 
field. The real vortices decay rapidly under the influence of viscosity (no vortex in the 
far field survives more than a few cycles). The purpose of the numerical calculations 
was not to simulate the effects of viscositv but rather to show that the behavior of flow 
(e.g., the evolution of the transverse vortex street) is primarily an inviscid phenomenon, 
dependent on the relative amplitude of the flow or the Keulegan-Carpenter number. 
The results show that the effect of viscosity is important in bringing about the 
separation and in decaying the vortices but not in bringing about the particular vortical 
structure. 

The vortex motion about the cylinder is better illustrated through the use of 
the vector plots of the velocity field. Figures 1.8 through 1.20 show the velocity field 
at suitable time intervals over one and one-half cycles, starting at T* = 0.5. At the 
end of the first half cycle (Fig. 1.8), two large vortices are separated by a strong current 
between them. As the flow reverses (Fig. 1.9) and the ambient velocitv reaches a value 
of |U| 


moves to the left. partly over the top and partly over the bottom of the cvliideragiiiae 


0.707, the vortex adjacent to the cylinder splits into two unequal parts as it 


vortex to the right of the current in Fig. 1.8 1s simiply convected towards the left of the 
cvlinder. At T* = 0.75 ({L| = 1), the vortices convected over the top of the eylmimes 
give rise to earlier separation on the upper left side of the cylinder. Furthermonemae 
lew vortex growing on the lower left side of the cylinder is pulled towards the top of 
the cvlinder. Figures 1.11 (JU| = 0.707) and 1.12 ({L] = 0) show the development of 
two large vortices, again separated by a strong current (compare Figs. 1.8 and 1.12). 
Figure 1.13 ({L| = 0.707) shows that the vortex adjacent to the left side of the ¢y iit 
remains nearly intact and is convected towards right over the top of the cvlinder. 
Figures 1.14 (T* = 1.25, JU] = 1) and 1.14 (T* = 1.375, |U| = 0.707), shew tigi 
separation of flow on the right side of the cylinder gives rise to new vortices, in a 
manner similar to that shown in Fig. 1.10. Subsequently, this leads to two large 
vortices (Fig. 1.16), separated by a current as in Fig. 1.8. Evidently, Figs. 1.3 ance 
are not expected to be identical in an ever evolving unsteady flow. Nevertheless, the 
essential features of the two flows (separated bv a full cycle) are quite comparable. 
The secondary vortical structures seen in Fig. 1.16 are a consequence of the inviscid 
nature of the calculations. As noted earlier, they do not survive the effects of viscosity 
and turbulence in physical experiments. Figures !.17 through 1.20 show the evolution 
of the flow for an additional half cvcle. Clearly, the flow pattern repeats itself with 


imnor differences (compare Figs. 1.8 and 1.16 and Figs. 1.12 and 1.20). 


It is evident from the foregoing that the reason for the occurrence of the 
transverse vortex street in the range 8<K<13 (corresponding to the amplitude-to- 
diameter ratios of A/D = 1.27 to 2) is, as expected, the formation of two asvmmetric 
vortices during a given half cycle and the effect of the returning vortices on the 
formation of the new vortices. 

The vortical structures and the transverse current shown in Figs. 1.8 through 
1.20 cannot be directly compared at an arbitrary time T* with those shown in Fig. 1.1 
since the streaklines of an oscillating flow about a stationarv cylinder are not the same 
aoetnose of an oscillating cylinder in a fluid otherwise at rest. Nevertheless, a 
Pommatison of Figs. 1.1D and !.1L with Figs. 1.12 and 1.20 {at times of zero ambient 
velocitv or zero cylinder velocity) has shown that the radial and the angular positions 
Siete vortices are predicted with surprising accuracy. In Figs. 1.1D and 1.1L the 
vortex to the left of the current is located at r/c = 5.45 and 9 (from the vertical) = 17 
degrees. The position of the same vortex in Figs. 1.12 and 1.20 is r/c = 5.42 and 0 
(from the x axis) = 15 degrees. It should be noted in passing that the side from which 
the domunant vortex sheds may become switched by the action of random disturbances 
in the flow or by stopping and restarting the flow. In the numerical calculations the 
direction of the transverse flow depends on the side at which the asymmetry 1s 
introduced into the flow. Had the asymmetrv been introduced on the lower side of the 
cylinder, the transverse flow would have occurred on the lower side of the cylinder. 

The creation, shedding and backward convection of the dominant vortex in a 
flow field, where there are only a few vortices, strongly affect the pressure distribution 
about the cylinder. Thus, it would have been desirable to compare the measured and 
predicted pressure distributions. However, no such data is available at K = 10 or at 
any other K value. Wilson (1983) measured the differential pressure between two 
diametrically opposed points on a cylinder and provided differential pressure 
distributions at suitable times. Calculations have been performed with the numerical 
model to obtain similar differential pressure data. Figure 1.21 shows two comparisons 
of the measured and calculated differential pressure distributions. Even though the 
magnitudes differ somewhat, the general shape of the two Cp curves are remarkably 
similar. The reason for the difference in their magnitude is directly related to the 
strength of the vortices. Previous investigations (see e.g., Sarpkava 1973) have shown 
that the calculated vortex strengths are somewhat larger than those obtained 


experimentally. An additional loss of vorticity could have been introduced into the 
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niodel in a manner simular to that used by Sarpkaya and Shoaff (1979) in order to bring 
the strengths of the shed vortices, and thereby the differential pressure distributions, 
into closer agreement with those obtained experimentally. However, it was decided to 
avoid a relatively subjective dissipation mechanism. It appears that the kinematics of 
the flow does not strongly depend on the strength of the shed vortices. This is in 
confornuty with the previous applications of the discrete vortex model (Sarpkaya & 
Shoaff 1979). 

Figures 1.22 and 1.235 show the calculated drag and lift coefficients as a 
function of time. The effect of the passage of a vortex over the Shouldepegigam. 
cylinder 1s exhibited in these figures in two ways. First, it brings sharp changes in Cd 
near its maximum. Second, it gives rise to a net transverse force on the cylinder. 
Evidently, the lift force 1s very sensitive to the flow asymmetry and depends strongly on 
the proximity to the cylinder of the vortex passing over the cylinder. The calculations 
wil have to be carried out over many more cycles in order to obtain a root-mean- 
square value of the lift coefficient for comparison with that obtained experimentally. 
This has not been done in the present study for the reasons cited earlier. 

A comparison of the measured and calculated drag force is shown in Fig. 1.24 
The agreement is not as good as expected primarily due to the fact that the strength of 
the calculated vortices 1s somewhat larger than those encountered in the experiments. 
A 5 to 10 percent decrease in circulation could have brought the measured and 
calculated forces into much closer agreement. As noted earlier, however, this was not 
the purpose of the present investigation. 

The attention will now be directed to the motion of the stagnation and 
separation points and to the vorticity fed into the shear layers. Figures 1.25 and 1.26 
show the motion of the upstream (primary) and downstream (secondary) stagnation 
points, respectively. Here upstream mieans the side of the cvlinder facing the ambient 
low. The primary and secondary stagnation points are approximately 1SO degrees 
apart, as would be expected. For example, at T* = 0.5 the primary stagnation point 
jumps to nearly zero and the secondary stagnation point jumps to approximately 180 
degrees. The spurious oscillations in Figs. 1.25 and 1.26, particularly in the position of 
the secondary stagnation point, are due to the disproportionately large influence of the 
proximity to the cylinder of a few point vortices, 

Figures 1.27 and 1.28 show the motion of the primary and secondary 


separation points (see also Fig. 1.3). Aside from the spurious secondary oscillations 


(due to the reasons cited earlier), Fig. 1.27 shows in general that the primary separation 
angle increases rapidly as the velocity of the ambient flow increases from zero (e.g., at 
ae —wiece 12 1) to aneaverdee value or about I10 degrees at [* = 1.5 (sce Fig. 
Into Ilen itedecreases rapidly to zero as T* approaches 2 (see Fig. 1.20). A 
comparison of Figs. 1.27 and 1.28 shows that there is a phase shift of one half of a 
cvcle between the positions of the primary and secondary separation points, as would 
be expected. There is, at present, no separation point data for comparison with those 
obtained experimentally. 

Finally, Fig. 1.29 shows the vorticity fed into each shear layer. The first and 
second shear lavers have relatively large and oppositely-signed vorticity. They are not 
each other's mirror inmiage due to the asymmetry of the flow. The vorticity in the third 
and fourth shear lavers is about 25 to 30 percent of that in the primary shear layers. 
Clearly. the vorticity fed into the primary shear layers (particularly that fed into the 
first shear laver) has a sinusoidal character to it. This is because the vorticity is related 
to the ambient flow velocity (sin@t) through the velocity UL, at the separation point. 
The flow variables such as vorticity cannot be measured directly. Thus, the accuracy 
of its magnitude and time-dependence can only be inferred indirectly through the 
comparison of the measured and calcuiated forces and pressure distributions. It 
appears, on the basis of such comparisons that, the calculated vorticitv is about 10 
percent larger than that prevailing in physical experiments. The reasons for this are 
partly the diffusion of vorticity by viscosity and turbulence and partly the three- 


dimensional nature of the physical experiments. 


D. CONCLUDING REMARKS ON OSCILLATING FLOW ABOUT 
CYLINDERS 


The discrete vortex model used in the present investigation avoids many of the 
problems associated with the Eulerian finite difference and finite element methods but 
is subject to some of its own. Chief among these is the excessive computer time needed 
for the calculation of the convection of vortices and the difficulty to account for the 
effects of viscositv and turbulence. 

In the present simulation the discrete vortex model has been combined with the 
boundarv layer calculations and the positions of the separation and stagnation points 
have been calculated for a Keulegan-Carpenter number of K = 10. The results have 
accurately predicted the formation of a half Karman vortex street in the transverse 


direction. The calculated positions of the vortices were found to be in good agreement 
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with those obtained experimentally. The measured and calculated drag force and the 
differential pressure distributions showed reasonably good agreement. The measured 
in-line force was somewhat smaller primarily due to the fact that the vorticity fed into 
the shear layers has not been artificially-reduced to bring the measured and calculated 
forces into closer agreement. The results have also shown that the effect of the 
backward convection of a large vortex over one side of the cvlinder is indeed very 
pronounced on all the measured and computed characteristics of the flow. This is one 
of the most important reasons as to whv the Morison’s equation (see e.g., Sarpkava & 
Isaacson 1981) fails to represent the in-line force acting on the cylinder with reasonable 


accuracy, particularly in the range 8<K < 15. 
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Figure 1.8 Flow field at T*= 0.500. 
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Figure 1.9 Flow field at T*= 0.625. 
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Figure 1.10 Flow field at T* 
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Figure 1.11 Flow field ar == 07375; 
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Figure 1.12 Flow field at T*= 1.00. 
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Figure 1.13 Flow fieldvat f= —siize: 
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Figure 1.14 Flow field at T*= 1.250. 
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Figure 1.15 Flow field art.S= 57: 
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Figure 1.16 Flow field at T*= 1.500. 
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Figure 1.18 Flow field at T*= 1.750. 
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Figure 1.19 Flow field at 1*— ssc 
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Figure 1.20 Flow field at T*= 2.00. 
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Comparison of numerical and experimental differential pressures. 
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The calculated drag coefficient. 
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Figure 1.23 The calculated lift coefficient. 
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Figure 1.24 Measured and calculated drag force. 
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Figure 1.25 Upstream stagnation point. 
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Figure 1.26 Downstream stagnation point. 
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Figure 1.27 Primary separation points. 
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Figure 1.28 Secondary separation points. 
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Figure 1.29 Vorticity fed into the flow field. 
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If. DISCRETE VORTEX ANALYSIS OF UNSTEADY FLOW 
ABOUT CAMBERED PLATES 


A. INTRODUCTION 

The determination of the deployment sequence of an axisymmetric porous 
parachute and the unsteady aerodynamic loads acting on it present a very complex 
coupled problem. The development of an analytical or numerical model which takes 
into account the effects of porosity, gaps, and variable opening schemes would allow 
numerical experiments on a large class of parachutes, reduce the number of expensive 
field tests to a few judiciously selected ones, and enable the designer to calculate the 
time history of the fall of the parachute and the strength required to survive the 
aerodynamic loads. However, the development of such a model is hampered bv a 
number of difficulties. 

The previous models for parachute loads are based by and large on empirical 
assumptions {see e.g., Heinrich and Saan 1987; Mcwey 1972). They rely on the 
observation that families of parachutes open in a characteristic length and seem to 
have aerodynamic properties that relate well to the projected area of the parachute. 
The apparent mass is assumed to be a function of the projected area only and is not a 
function of the prevailing flow characteristics. The vortex sheet analysis was used by 
Klimas (1977) to derive the acceleration-independent apparent mass coefficient for 
arbitrary-shaped axisymmetric surfaces. Muramoto and Garrard (1984) used a 
continuous-source model to predict the steady-state drag of ribbon parachutes. The 
analyses did not, however, deal with the evolution of the unsteady wake and its 
interaction with the canopy. 

Itis in view of the foregoing that a fundamental study of the separated time- 
dependent flow about two-dimensional rigid cambered plates was undertaken. Clearly, 
the flow about a rigid cambered plate is considerably simpler than that about a porous, 
axisymmetric, and flexible parachute and the results, regardless of the degree of their 
agreement with corresponding experiments, may not have direct relevance to the 
practical problem under consideration. But the object of this investigation was the 
understanding of the evolution of the wake under controlled conditions rather than to 
provide a design tool. It 1s hoped that an investigation of this type will reveal the 


underlving physics of the phenomenon (particularly that of the parachute collapse), 


help to interpret the full-scale results and will provide inspiration for the development 
of more general vortex models with which the dynanucs of axisymmetric, porous, and 
flexible parachute canopies can be investigated. Efforts directed towards the 
development of a general numerical model, driven by the ever-present pressures of 
practical considerations, are deemed somewhat premature. Such efforts will have to 
face not only the problem itself but also the deficiencies of the vortex modelsmame 
attempt to address to both of them simultaneously. The model presented herein 
removes the ambiguities associated with the use of the discrete vortex model and 


provides results which are in excellent agreement with those obtained experimentally. 
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IT. ANALYSIS 


A. TRANSFORMATIONS AND THE COMPLEX VELOCITY POTENTIAL 

The calculation of the velocity of any one of the vortices and the force acting on 
the body requires a conformal transformation (in which the camber becomes a circle), 
a complex-velocity potential representing the vortices, their images, and the two- 
dimensional irrotational flow around the body, and the use of the generalized Blasius 
theorem. 

The flow in the circle plane may be transformed to that about a cambered plate 
through the use of two successive transformations, one from ¢ plane to the €° plane 


and the other from the €° plan to the z plane. These are given by (see Fig. 3.1) 


a. 7 b. n° 
—> vane 6 a 
¢-plane SY ¢-plane 





Figure 3.1 Circle and physical planes. 
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z=CO°o-—, and €°= +m (3.1) 


Combining the two, one has a direct transformation from the ¢ plane to the z 


plane as 


> 


~~ 


7 = CP mee 





= (32) 
eer Til 


[t 1s easy to show that the camber in the z plane is a circular arc. 

The y-axis in the z plane passes through the tips of the camberaaamue 
advantageous to locate the origin of the coordinate axes at the geometric center of the 
camber, 1.e., at the center of the circle part of which represents the camber asus 


easily accomplished by shifting the origin of the coordinate axes by 


n= 1 
Z. = Bre) 


2 m 





Where Z, 1s the x coordinate of the origin of the circle in the z plan. Thus, one has 


-~ 


ao 





Z=CPm — se Ze NW Zh ates a 2 (3.4) 


eae ‘i 
which transforms the circle in Fig. 3.la to the physical plane in Fig. 3.!d. Table 1 
summarizes the relationship between m, Z;, the included angle of the camber, b, and 


the radius of the caniber. 


TABE Eel 
SUMMARY OF THE PARAMETRIC KEDAMIRG SS init 


m Zz. 2a b R=1/m 
cos60 = .5 -1 le 866 2 
Cc 55— ey 0. 180 0 /2 
cos30 =  S6Gumle 3 240 iS 2!./3 


The complex potential function W in the circle plane (see Fig. 3.la) which 
describes a uniform flow U (assumed to be time-dependent) with a doublet at the 
Origin to simulate the cylinder, Lia clockwise-rotating vortices (called q-vortices), I, 
counter-clockwise rotating vortices (called p-vortices), and the images of all the p-and 


q-vortices in the circle plane may be written as 


; co iT, iT, ¢? 
eG aa Ene = Sop) eet See) 
g 27 2 Op 
aD Ee ceil 
+ VY —KPpact - —~ VY —kepat - —) —- —9bwt - € 
2 n(§ — &.) er nf de nS — pq) 
iT, eel ue . 
Ponta — 2) — Seine - 4 + Feng - 35) 
q 
: nce Keen Skq 


in which es and Sp Tepresent respectively the strength and location of the k-th p- 
ORLEX, ie and Skq the strength and location of the k-th q-vortex, and c the radius of 
the cylinder; an overbar indicates a complex conjugate. The need for the separate 
identification of the p-and q-vortices and for the singling out of one of the vortices in 


each shear layer (namely Ty, and Dog: nascent vortices) will become apparent later. 


B. COMPLEX VELOCITIES OF VORTICES 

The convection of the vortices and the calculation of the forces acting on the 
body require the evaluation of the velocities at the vortex centers. For the velocities in 
ivermcitcle plane this reduces to subtracting from Eq. (3.5) the complex potential 
corresponding to the vortex for which the velocity poets. are COMmeedctcrimined 
and evaluating the derivative of the remaining terms at § = ¢,. To determine the 
velocities in the physical plane, however, one has to subtract bie fe En Z ee irom 


ie().5) OF, i terms of € , the terms (see e.g., Sarpkaya 1967, 1975) 


bh 








Cn (qae en lier (3.6) 


Ln|1 + —__ 
(¢+m)(6, +m) 
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It should be noted that the first term in Eq. (3.6) ts the complex finenen 
corresponding to the k-th vortex in the § plane. The second term appears merely as a 
consequence of the transformation used. 

The above procedure may be generalized as follows. Consider the potential 
function for a single vortex in the physical plane and ignore, for the time being, the 


multiplier in front of the logarithmic term (ie., if,/2m). Then one has 


Ln(z — z,) = Lo{ {6) — €6,)]  withz = 0) (3.7) 


equation (3.7) may be written as 


(O) — AS) 


LIn(z-z)=In(G-G6)+L 
n(~ 2) = In ~ &) + La 


(3.8) 


Evidently, the first term on the right hand side of Eq. (3.8) represents the vortex in the 
circle plane. Let us now examine the derivative of the second term with respect to z. 
One has, 
=] 
dw, ta /aO-Ko)\/no-as) \ ae aa 
dz dc oe Bornes dz 


where dz'd¢ = f(t). In the neighborhood of ¢. the function f (¢) may be expanded as, 


RG) Gerona ae f(a (3.10) 
Thus, one has 

dW | 1 aaa oa 

—F= Lim ——__—_—_—_—_— |-— . 

dz Ce i) ST) 0S) 


fac 


Or 
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dW Fo.) 





ap So (a2) 
dz 2F* (G,) 
Thus, the complex velocity in the physical plane reduces to 
; 4 
Pee tate) |- —_ (3.13) 
=U iv =— —- —Ln(o - —-|- — = Bs 
x 7 de : 2n 2 OW FQ} 4 £206.) 
oe 5, 
in which for a p-vortex 
b? 
fre) (3.14) 
e (cc m) 
and 
; ope 
fgj=-— (als) 
Mgeniast term in Eq. (3.13) reduces to 
2 
- Le Ga Morn 16.0 (3.16) 


2x [G+ m)- + b*? PF 


Mase result could have been deduced directly from Eq. (3.6). However, the 
generalization of the method enables one to apply Eq. (3.13) to any vortex for any 


transformation between the circle and the physical plane. 


C. KUTTA CONDITION 
The fact that the flow separates tangentially with a finite velocity at the edges of 


the plate (Kutta condition) may be expressed by requiring 


— =0 SS i Greene 1 (eel) 
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Thus, insétting Eq. (3-5)n eq, ay emnesna. 
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cL 2 t ¥ 
SKp Sk 
Cc? 
—U{ l= == 0. (3.18) 
or 
Equation (3.18) may be decomposed into two parts as 
Ps iP ae I = To, a l 
w 9 w 
2m b. ~ Sop C = 20 G, ~ $0q i ac 
0p . S0q 
tel Se ee (3h 


where the terms containing the strength of the nascent vortices represent the velocity 
induced at the tip of the camber by the nascent vortices and the term in parenthesis the 
velocity at the tip due to all other vortices (and their images), the doublet at the center 
of the circle in the ¢ plane and the ambient velocity. 

Equation (3.19) represents two coupled equations for the strengths and positions 
of the nascent vortices. Thus, the solution of the said quantities does, in general, 


require an iteration. However, this iteration mav be avoided by noting that the 


velocity induced by a nascent vortex at the opposite tip is very small and certainly 


negligible. Thus, Eq. (3.19) for one of the nascent vortices may be reduced to 








IF, l l Z 
—- —" —_—- Ss yt (et vee (3.20) 
2m a7 Soq - oe 
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A similar expression may be written for the other nascent vortex. The use of the Kutta 
condition, as expressed by Eq. (3.20), will be further explained following the discussion 
of the tip velocity. It suffices to note that all nascent vortices satisfying the Kutta 
condition do not yield either the same tip velocity or the same velocity distribution in 
the neighborhood of the tip. There are, in fact, certain preferred positions for the 
nascent vortices which yield physically realistic velocity distributions near the tips of 


the cambered plate. These nascent vortex positions will be discussed later. 


D. TIP VELOCITY 

According to the Kutta condition the tangential velocity at the tip is finite. The 
purpose of the following is to determine this finite velocitv. It may be determined 
either through the use of l’Hopital’s rule or through the use of a more general 
expression which is valid for all other transformations. 


The velocity at the tip is given by 


dw dw de 
— = — + at Z 
dz dt dz 


zx 2ib (3.21) 


For an arbitrary point z, Eq. (3.4) yields, 


dq l S10 
— = — + ————_— (3.22) 
dz Z heen Z, | 
In general, one may write Eq. (3.22) as 
d I (nae? 
dz 2 A (emmziy= t AD= 
or multiplying both sides with ./(z — z,), one has 
dg 
a= (2-2) * = 
ie (ZZ) i 
l ane 
= —(z —z)'* + eee (3.24) 


2 aa eb) 
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Thus, as Zn = 216, one has 


os See a 
f(z) = (z-2z)'% = = mae (3.25 
or 
dc Jib Die 
= == = v (3.267 
dz aV(Z57Z,) Gomes, 
Expanding dW/d in the neighborhood of ¢ one has 
dW ; Ww 
= (GCN) ae (3.27) 
Combining Eqs. (3.26) and (3.27) one finally has 
(3.28) 


dw d-W , 
ae ea 2) 


pe = oo t 
dz ge 


t 


= ee : ‘ 2) are: 
Noung that for the case under consideration f (z,) = 16/4, one has 


dw] d°W ae, a 
—_—|= —s- (— J.2 
dz d¢- 2 





t 
Equation (3.29) yields the desired finite tip velocities. It 1s easy to show that it may be 


obtained directly from Eq. (3.21) through the use of l’'Hopital’s rule. 


E. TIME DEPENDENT-FORGES 
The force acting on the body in the physical plane may be calculated either 


through the use of the pressure distribution or through the use of the rate of change of 


impulse. 
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Bernoulli's equation for unsteady flow 1s given by 


as Py et oo ee 
—+—1)-(— + —) - f —4s = 30 
or are cae 7 (3.30) 


where the indices indicate two points on the body in the physical plane. Since there ts 
no pressure drop across the shear layer and since the integral term in Eq. (3.30) 1s zero 


at the tip (1.e., ds=0), one has 


Nees eV” 
fn) = a ~ a (3.31) 


where V,, and V,, represent the tangental velocities on the upstream and downstream 
faces of the tip. It 1s important to note that f (t) in Eq. (3.31) is also the time rate of 
change of circulation, 1.e., the rate at which vorticity is shed into the wake from the tip 
of the cambered plate. 

The normalized form of Bernoulli's equation between any two points m and n 


then becomes 


cP Vite ame Vo 
Ht = tL te + ee a7 Ceo 
pe; 2 U U 


oO Oo 


The integration of the differential pressure between the upstream and downstream faces 
of the camber yields the force components in the x and y directions, i.e., the drag and 
lift forces. 

The force acting on the body can also be calculated through the rate of change of 


impulse. It is given by 


Z 


= 210 ae g 2 
= 4npc*U (1 — ey + a r Za a Zi) | (3:35) 
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which may be written as 








C.+iC¢-—+= ie 
= * = N - 2 ae = ae 
dL 29U 2b Ube Dee 
Cr We i é c 
+ — ——_-_——_ {[_—_ Sf (=) 3334 
% AUA0 Te (= = (3.34) 


in which U, is the reference velocity, U, the rate of deceleration of flow and Zee 
1.e., the transformation given by Eq. (3.4). Equation (3.34) may also be deduced 
directly from the generalized Blasius equation. It is important to note that the force 
calculated from Eq. (3.54) includes the effect of the rate of change of circulation 
between two successive time steps. Thus, it may be smaller or larger (depending on the 
sign of I) than the force calculated through the integration of the instantaneous 
differential pressure Eq. (3.32). This is because of the fact that the instantaneous 
pressure depends only on the prevailing flow conditions and does not account for the 
rate of change of total circulation between successive time steps. In the calculations to 


follow L. and c are taken as unity for sake of simplicity. 


FL. METHOD OF CALCULATION 

The methods used in the past in the determination of the vorticity Mux@iven 
sharp-edged bodies may be roughly classified into two broad categories. The first of 
these involves the use of variable nascent vortex positions (see e.g., Sarpkaya 1968, 
1975) and the second, the use of fixed nascent vortex positions (see e.g., Clements 
197521975): 

The method of fixed positions involves the selection of a suitable fixed point in 
the flow near the separation point and the use of the velocity U, at that point to 


calculate the rate at which vorticity is shed into the wake from 
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In this method the positions of the nascent vortices are the crucial parameters. I[he 
previous applications of this method did not examine the effect of the position of the 


nascent vortices on the velocity distribution in the neighborhood of the separation 
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point. Only the distance of the fixed point to the body was varied and bracketed 
between two subjective limits by comparing the calculated results with those obtained 
experimentally. In this method no interaction is allowed between the shed vortices and 
the amplitude of oscillation of the point or the time of appearance of the nascent 
vortices. Furthermore, the time interval is chosen more or less arbitrarily (Kiva and 
Arie 1977) (repeating a few calculations with a single program with only the time step 
changed and also by referring to the results of the previous investigations). Thus, the 
velocities at the outer edges of the shear layers are only indirectly related to the 
strength of the nascent vortices and the fixed time interval. Evidently, the velocities in 
[me immer and outer edges of the shear layers, the time interval, the strength and 
position of the nascent vortices, and the Kutta condition are interdependent and that 
both the position of the nascent vortices and the time interval cannot be chosen 
arbitrarily, even if they are chosen judiciously on the basis of previous experience and 
trial calculations. 

Sarpkaya (1975) used the method of variable nascent vortex positions and 


deternuned the rate of shedding of vorticity from the relation 


roo. 
a 


BoC 
Where U., is interpreted as the velocity in the shear layers calculated by using the 
average of the transport velocities of the first four vortices in each shear layer. The 
positions of the nascent vortices are chosen so as to satisfy the Kutta condition at the 
edges of the body and thus they can move slightly with time. Thus, this method 
simulates in a satisfactory manner the mechanism of feedback from wake fluctuations 
to the fluctuations in the rate of circulation. The number of disposable parameters is 
reduced to a minimum and in this sense this method is superior to the method of fixed 
positions. However, the use of the average of the transport velocities of the first four 
vortices remained questionable. 

It was often assumed that the vorticity flux could not be calculated, at each time 
interval, as it 1s applied to sharp-edged bodies, through the use of the mathematically 
finite velocity occurring at the sharp edges of the body. This assumption was based on 
the fact that the separation points are singularities of the transformation used and the 


numerical procedures may not be stable. 
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It is on the basis of the foregoing that an original study was undertaken to 
establish once and for all a method whereby the nascent vortices may be introduced 
into the flow without any ambiguities. The method finally arrived at will be explained 
through the use of a series of figures and velocity plots. 

Figures 3.2a and 3.2b show the tip region in the circle and physical planes, 
respectively. The regions A and B in Fig. 3.2a were discretized through the use of a 
suitable grid and a single vortex was placed at a grid point. The strength of the vortex 
was determined from the Kutta condition Eq. (3.20). Then the velocity normal to the 
radial line OZ in the physical plane (Fig. 3.2b) was calculated in the vicinity of the tip 


through the use of the complex velocity potential. 





Figure 3.2 Tip region in the circle and physical planes. 


Placing the vortex along the radial line OM (in the € plane) vields a single valued 
tip velocity independent of the strength and the position of the vortex and dependent 
only on the plate geometry, 1.e., b and 0. It is easy to show that the velocity at the 
edge of the plate reduces to q, = +(1b/2) e tis (— 1 —e “185) For the casemannd 
120-degree camber this gives an absolute value of 0.433 with a velocity direction 


s 


Opposite to that expected at the edge of the plate (siigmaca) 
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Figure 3.3 Velocity profile along the radial line (nascent vortex on OM), 
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Placing the nascent vortex to the right of the radial line OM in the circle plane 
(or along the circular arc in the physical plane) always requires a stronger vortex to 
satisfv the Kutta condition and results in a tip velocity which is unrealistic both in 
magnitude and direction (Fig. 3.4). j 

Placing the nascent vortex along the radial line OZ in the physical plane or 
outside the region A shifts the point of maximum velocity away from the edge of the 
plate (towards the downstream side). This, in turn, results in a leakage of fluid through 
the shear laver and requires a stronger vortex to satisfy the Kutta condition (Fig. 3.5). 

The entire region A enclosed by the transformation of the radial line OZ in the 
physical plane and the radial line OM in the ¢ plane, is examined to determine the 
most appropriate positions of the nascent vortex. Figures 3.6 through 3.8 show a 
three-dimensional plot and the contour lines of [, U(max),U(tip), and U(tip) as a 
function of the radial positions R and the angular positions R@ for those locations of 
the nascent vortex for which F <1, U(max)'U(tip)=<6, and Up) et 

Figures 3.9 through 3.11 show three representative velocity profiles fomermmes 
different positions of the nascent vortices in the region defined above. The most 
striking feature of these figures is that the maximum Velocity near the tip can exceed 
considerably and unrealistically the velocity at the tip and that only for certain vortex 
positions does the maximum Velocity (the velocity on the inner face of the camber at 
the tip) approaches smoothly the finite tip velocity. These calculations have shown 
that there 1s, in fact, a finite region in which the nascent vortices may be introduced in 
order to produce a tip velocity which is nearly equal to the maximum velocity in the 
vicinity of the tip. Clearly, it is only for unique combinations of the radial distance R 


and the angular position R@ that the said velocity ratio is equal to unity. 
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Figure 3.4 Velocity profile along the radial line (nascent vortex to the right of OM). 
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Figure 3.5 Velocity profile along the radial line (nascent vortex on OZ), 
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Figure 3.6 Three-dimensional plot of and the contour lines for I distribution. 
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Figure 3.7 Three-dimensional plot of the velocity ratio. 
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Figure 3.8 Three-dimensional plot of the tip velocity. 


719 


° ° 
e 0 ° 2 
Oo oO Oo ° 
° 9° ° ° 
° x 2009, ° 56 = 
0900,.9 0° 950 ° 
oR 0 %9 |° 
5 
/ 
/ 
/ 
/ 
4 
‘ GAMA - -0.70 
ae 
e 
@ @ 
e 
Ww —< 
- e 
N se" 
e 
SG 





0.875 0.315 ‘ tans 


Reine 
“ 
me eal 


Figure 3.9 Velocity profile along the radial line (nascent vortex at €=0.05, 60 = 25.0°). 
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Figure 3.10 Velocity profile along the radial line (nascent vortex at €=0.10, 68 = 25.0°). 
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Figure 3.11 Velocity profile along the radial line (nascent vortex at €= 0.20, 60 = 25.0°). 


Figure 3.6 shows the required vortex strength as a function of a R and RO. Even 
though these calculations have been carried out with a single vortex, the subsequent 
calculations with larger number of vortices have shown that the relationship between 
the maximum velocity and the tip velocity remains practically invariant as long as the 
nascent vortices are alwavs introduced at the fixed point which produces 
inde (tipje = sienor the sinele vortex. The reason for this 1s that the said velocity 
ratio is primarily dictated by the nascent vortex and the complex potential used in Eq. 
(3,29). to calculate the tip velocity, encompasses the effect of all the vortices in the flow 
field. In other words, it does not make any difference whether the velocity q..1.e.. (U, 
+ iv.) in Eq. (3.20) is produced by the ambient velocity or by a large number of 
vortices in the field. 

The fact emerging from the foregoing analvsis 1s that the nascent vortex cannot 
be placed arbitrarily (e.g., along the radial line in the circle plane or along the extension 
of the camber in the physical plane). Otherwise, the velocity distribution in the vicinity 
of the tip becomes unrealistic and unrepresentative of the evolution of the shear lavers 
emecmner side Of the camber. Furthermore, one is then forced to make arbitrary 
assumptions regarding the strength and the convection of the nascent vortices. 

The foregoing extensive analysis led to the conclusion that the nascent vortices 
should be introduced at r = 1.0925 and 9=0. + 2.08°, for the case of the 120-degree 
camber. To be more precise, the two nascent vortices are placed at the angular 
positions 0, Uo) weaves rand 0, = 0 ,+ 2.08, during the period for which V, - V, > 
oe or Y, — V,, the positions of the two nascent vortices are switched to their 
Sonresponding images with respect to the radial line OM, 1.e., they are placed at 0 = 
Q., + 2.08 and 0,3 = §.,— 2.08. The evolution of the very early stages of the flow in 
ihe immediate vicinity of the tips of the camber is shown in Fig. 3.12. 

Pio edetine cloctty field 1s a consequence of the first two nascent vortices 
Maeroduced at the points noted above. Figures 3.12b through 3.l2d show the 
development of the flow field and the starting tip vortex subsequent to the introduction 
of the 4%, 8", and 11" nascent vortex, respectively. 

There 1s not a unique procedure for relating the rate at which vorticity 1s shed 
into the wake, the Kutta condition, the velocity with which the nascent vortices are 
convected, and the time interval for the convection, all of which help to simulate the 
Experimentally observed features of the free shear layers. Fage and Johansen (1928). 


through quite ingenious experiments with steady flow about various bluff bodies, have 
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shown that vorticity is shed from the two sides of an axisymmetric body (a circular 
cvlinder) or a sharp-edged body (a plate normal to the flow) at the same rate: that the 
motion in a sheet is steadv near the body, except possibly near the inner edge of the 
shear layers; that fluid flows into a sheet through both edges, but at a greater rate 
through the outer edge; that at each section of the sheet the velocity rises from a small 
value to a Well-marked maximum value (approximately Se eee) amid then: ver 
slowly decreases to about 1.35 within a distance of approximately v~ 2c, where the 
breadth of the sheet reaches a value of A ~ c; and finally, that the velocity V, at the 
outer edge of the sheet is much larger than the velocity V, at the inner edge (except 
during the deceleration period of the flow) and ce may be ignored in calculating the 
vorticity flux from ¢T ct =0.5 ( Vie = VO 

Pimeiespresent calellawonetne “Ortex Sttemeih, the velocities On either side of the 


Shear laver, and the time interval are related by 


ore | 


lr = 0.5(V,° — V,)-At cen 
in Which I is the strength of the nascent vortex, V, = U(tip) and V, is the velocity at 
the downstream face of the camber near the tip. The velocity V, can be calculated 
Gorrectly in a numiber of ways, to be described later. Suffice it to note that in general 
V. 1s very small (for steady flow) and that the method of its calculation has very little 
Seno iniiiience on the strength of the nascent vortex or on the time interval to be used 
for a given vortex strength. The velocity V, becomes important only when the wake 
begins to move towards the camber (1.e.. during the period of flow deceleration). 

To explain the computational details of the method let us consider a particular 
time t after the start of the motion and assume t to be sufficiently large so that there 
are a number of vortices in the wake. Then the appearance and convection of the 
vortices proceed as follows: 

oe Determine the strength of the mascent vortices from the Kutta condition Eq. 
(6520) ine niet So is a known fixed position for each nascent Vortex ; 

mie Place the nascent vortices at Sop and S04 and calculate the velocity V, at the 
two edges of the plate ; 

(3) Calculate V,, representing the velocity at the inner boundary of the shear 
laver, as the average of the velocities at three points along the radial line OZ 


in the physical plane i.e., atr = 1, 1.05, and 1.1; 


(4) Calculate the ume interval, for eacn edge of the plate, from Eq. (o.3.))ememm 
tne known values of [ and the velocities V, and V,. Store the average of the 


two time intervals for use in subsequent calculations; 


oy 


(5) Calculate the velocity induced at the center of all other vortices; 

(6) Convect the two nascent voruces with a velocity 0.3(V, + V, ) for an average 
time interval At (note that the vorticity 1s convected with the averacewveremnes 
of the shear laver). If the distance travelled by a nascent Vortex 1S nomena 
0).03=0.01, 1 1s convected twice for a time interval At 2. The Sudseqmiems 


convections of the nascent vortices are made using the velocity ineWeequam 





(heir center 
(7) Convect all other vortices for the same time-interval At using a second order | 


scheme given by 


AtrAt) = z(t) + 0.5 [3zt) = "Zit Ae ( 


Gd 
ti) 
CO 
areas 


(ane lee — tee ie 

(S) Remove the vortices from the calculation whenever they come nearereumem 
Q).03 to the camber in the physical plane (except the first 20 vortices Irommemmm 
olen 

(9) Coalesce the same sign vortices with a separation of less than 0.03 (inwiae 
pnvsical plane. except the first 20 vorices): 

({0)  Catculate the tangentia! velocities and pressures on the inner and@vourermaees 
of the camber. Determine the drag and lift forces througn the intégmacieraes 
pressure and through the use of the rate of change of impulse. ake plotsven 
suitable variables (e.g., velocity distribution near the tip, variation of nascent 
voriex circulation with time, evolution of the wake, etc.): 

(11) Check the flow conditions to determine the state of the calculations: 

(a) IfV,—V, > 0.2 repeat the foregoing steps, 

(b) Stop the introduction of nascent vortices if 0 < V, = VG seieune 
RE tC ees Ou NOeuey 

(c) IfV, > V, switch the angular positions of the nascent vortices to their 
image points. Calculate V, as the average of the three velocities, at the 
upstream side of the tip of the camber, at three radial locations (0.93. 


0.9, and 0.83) and repeat the roregoing stemoaname 


S6 


(12) Make plots of the variations of various flow parameters (e.g., tip velocity, 
nascent vortex circulation, evolution of the wake, force coefficients, etc.) and 
terminate the run. 

The foregoing steps are quite general and can be used for anv camber, provided 
that the optimum points of placement of the nascent vortices are determined through a 


similar analvsis for the desired camber angle. 
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[V. DISCUSSION OFRESULRS 


A. NUMERICAL AND PHYSICAL EXPERIMENTS 
The calculations Were carried out for a time-dependent normalized velocity given 
by 


tt 
for T*= 
E 


x | . 
iI 
1A 
i 
ON 
Ca 


(4.1) 


and 


lI 


L— 0.1539(T* —8.65)+ 0.00531(T* — 8.65)? (1 


Oo 


which corresponds to that encountered in a series Of eXperiments carried (OU 
Vertical water tunnel. A detailed description of the equipment and procedures 1s given 
by Sarpkaya and Ihrig (1986) and will not be repeated here. Evidently, the caleulaiigms 
can be carried out for any specified variation of the velocity. For the casemungen 
consideration, the flow begins to decelerate at T* = 8.65 and the velocity omumme 
ambient flow reduces to zero at about T* = 19 (seeumiceea 

The computer program provided. at times specified, the positions of all the 
vortices, the rate of shedding of vorticity from the tips of the camber, the vyeloe@ins 
distribution on the upstream and downstream faces of the camber, the total and 
differential pressure distributions, and the force coefficients. 

Figures 4.2 through 4.4 show, at T* = Ujt'c = 4.35, the evolution of the wake. 
the tangential velocities at the upstream and downstream faces of the camber, the 
velocity profile along the radial line passing through the tp, and the total and 
differential pressure distributions (these plots are available at every time step but are 
not reproduced here for sake of brevity ). These and other figures show that the 
characteristics of the flow develop symmetrically prior to the onset of decelerauon (I* 
< §.65) and the differential pressure is positive everywhere (1.e., the pressure Insideumme 


camber 1s larger than that outside). 
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Figure 4.1 Variations of the velocity and the acceleration. 
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Position of vortices and tangential velocity distribution. 


Figure 4.2 


90 





S/O Sse .Se 


I 


1 


Sc) | 


OS‘ I 


Gob -9/1N = 


T 


ova 


ZY SNOW 31WsSOY¥d ALIDONA 


Figure 4.3 Velocity profile along the radial line passing through the tip. 
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Figure 4.4 Pressure distribution on the upstream and downstream faces. 


Following the onset of deceleration (see e.g., Figs. 4.5 through 4.7 at T* = 
10.84). the differential pressure near the axis of the camber becomes increasingly 
negative. The reason for this is that the deceleration of the flow brings the vortices 
closer to the camber. The significance of this result is that had the model been flexible 
(as in the case of a parachute) the central part of the camber would have collapsed as a 
result of the particular deceieration it is subjected to. Evidently the collapse 
phenomenon would not have remained symmetrical, as evidenced by field experiments 
with axisymmetric flexible parachute. Furthermore, it would have required the analysis 
of the flow about a flexible camber. It is because of this reason that in the present 
analvsis the camber is assumed to remain rigid. 

Pore SlATcermraaeavoltr lS ({Or the ambient flow Under consideration), the 
velocities induced at the downstream edges of the camber by the large vortices moving 
sidewavs and towards the camber give rise to oppositelv-signed vorticity. This, in turn, 
leads to the rapid growth of the secondary vortices (see e.g., Figs. 4.8 through 4.10 at 
T* = 17.86). The secondary vortices are relatively weaker than the primary vortices 
partly because thev have been in existence only for a short time and partly because the 
vorticity flux is not as large as that in the primary shear lavers. Consequently , the 
Semenotd Of the secondary vortices tends to orbit about the centroid of the primarv 
mOrtices. 

A comparison of Figs. 4.7 and 4.10 shows that the region of negative differential 
pressure grows with time and occupies a large central portion of the camber. In fact, 
the drag force acting on the camber becomes negative, as it will be seen shortly. 

Ficires 411 through 4/16 show, at suitable times, the velocity ficld about the 
camber. The rapid growth of the wake during the period of steady uniform flow 1s 
exhibited in Figs. 4.11 and 4.12. Figure 4.13 nearly corresponds to the time at which 
the deceleration is imposed on the flow. Figures 4.14 through 4.16 show clearly the 
backward motion of the primary vortices and the rapid growth of the secondary 
vortices. It is seen from Fig. 4.16 that the fluid motion is entirely due to the motion of 
the vortices in the flow field. The two vortices on each side of the camber form a 
counter-rotating couple and remove themselves rapidly from the field under the 
influence of their mutual induction velocity. Subsequently, the absolute value of the 
the differential pressure begins to decrease. Eventually, the differential pressure reduces 
to zero everywhere on the camber as the conditions approach to that of a bodv ina 


fluid at rest. 
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Figure 4.6 Tip velocity at T* = 10.84. 
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Figure 4.7. Pressure distribution at T* = 10.84. 
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Figure 4.8 Tangential velocity at T* = 17.86. 
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Figure 4.9 Tip velocity at I* 
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Figure 4.10 Pressure distribution at T* = 17.86. 
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The velocity field about the camber at I* 


Figure 4.11 
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Figure 4.12 The velocity field about the camber at T* = 6.05. 
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The remainder of the discussion of the results will deal with the variation with 
time of the tip velocities and the drag and lift forces. 

Figures 4.17 and 4.18 show the velocities V, and V, as a function of T*. The tip 
velocity V, decreases from an initially large value of about 3.5 to a nearly constant 
value of about 1.5 just prior to the onset of deceleration. Subsequently, V, decreases 
rapidiv during the period of deceleration and prior to the inception of the secondary 
separauon. Then V, increases to about 2 because of the backward motion of the large 
vortices near the tps of the camber. Finally, V, decreases once again as the primary 
and secondary vortices move sideways and away from the tips of the camber due to 
their mutual induction (see Fig. 4.16). 

The variation of V, with T* is significant only during two, relatively short. ume 
Intervals: at the start of the motion and at the start of the deceleration. These are the 
periods during which the vorticity flux changes rapidly in order to maintain the Kutta 
condition. During the remainder of time V, is negligibly small, as expected on the 
basis of the pioneering experiments of Fage and Johansen (1928) with steady flow over 
Various tvpes of bluff bodies. 

Figures 4.19 and 4.20 show the variation of the drag and lift coefficients as a 
function of time. The former is based on the integration of pressure and the latter on 
the rate of change of impulse. The drag coefficient calculated through the use of the 
rate of change of impulse is somewhat larger than that obtained through the 
integration of the instantaneous differential-pressure distribution. This is due to the 
fact that the impulse expression includes the rate of change of circulation between 
successive time steps whereas the pressure expression does not. I[t is a well-known fact 
that in real fluids the memorv of the fluid resides in its vorticity. Whereas in inviscid 
flows there is no memory and the dynamic characteristics of the flow (pressures and 
forces) are functions of only the instantaneous state of the flow. The adinaieam 
presented herein is for an inviscid fluid even though the phenomenon concerns the 
motion of a real fluid. The question of whether the rate of change of circulation 
should be included or excluded in the discrete vortex analvsis (first discussegmis 
Sarpkaya in 1968) is an unsettled issue. It appears that only the comparisons with 
experiments can clarifV the question. 

Figures 4.19 and 4.20 also show that C, rises rapidly (due to the rapid 
accumulation of vorticitv in the growing vortices) and begins to decrease as the 


vortices develop under the influence of a constant ambient velocity. Then the force 
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Figure 4.14 The velocity field about the camber at T* = 14.20. 
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Figure 4.15 The velocity field about the camber at T* = 16.30. 
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.16 The velocity field about the camber at T* = 2] 
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17 Variation of V 
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Figure 4.18 Variation of V, with time. 
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decreases sharply at the onset of deceleration and goes through zero near the middle of 
the deceleration period (T* = 11). The force acquires its largest negative value 
towards the end of the deceleration period. Subsequently, the force graduaily decreases 


ie zero. . 

Also shown in Figs. 4.19 and 4.20 is the variation of the lift force. It is negligible 
even in the later stages of the motion. This is primarily due to the fact that there 1s 
not sufficient time for the development of alternate vortex shedding either during the 
period of steady flow or during the period of rapid deceleration. 

Figure 4.21 shows a comparison of the calculated (through pressure integration) 
and measured drag coefficients. In general the agreement between the calculated and 
Weermiecasured drag C@eliicient 1s quite good. In the time intervals between 13 and [6 
and between 19 and 22, the calculated C, 1s somewhat larger. The reason for this is as 
follows. In the said time intervals the drag coefficient is relatively small and the 
viscous effects are relatively important in dissipating the vortices. This is not taken 
into consideration in the numerical analvsis. It 1s possible to introduce a small 
artificial reduction in circulation in order to bring the calculated and measured values 
into closer agreement. This has been avoided in the present analvsis in order to keep 
the discrete vortex analysis as pure and simple as possible. Figure 4.21 also shows 
that the calculation of the drag coefficient through the integration of pressure is 
superior to that through the use of the rate of change of impulse. 

Finally. a comparison 1s made between the calculated and photographed flow 
fields at corresponding times. Figures 4.22 through 4.24 show at umes T* = 6.05, §.55 
and 16.30 the flow in the immediate vicinity of the camber (plotted to the same scale). 
shows that the agreement between the calculated and observed flow fields is indeed 


very good. 


B. CONCLUDING REMARKS 

The results presented in this section have shown that the discrete vortex model 
can be used with confidence to predict the evolution of the wake about a cambered 
plate immersed in a an arbitrary time-dependent flow. The evolution of the wake 1s 
remarkably similar (including the formation of the secondary vortices) to that obtained 
in flow visualization experiments. The drag coefficients resulting from the analysis and 
Eeperlmients agree reasonably well. Ilus agreement can be improved with the 
introduction of a small circulation dissipation. The drag coefficient calculated through 


the integration of the instantaneous pressure distribution agrees more closely with that 
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Figure 4.19 Drag coefficient calculated from pressure integration. 
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Figure 4.20 Drag coefficient calculated from impulse expression. 
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Figure 4.22 Comparison of flow fields at T* 
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Figure 4. 
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obtained experimentally and points out the fact that the inclusion of the rate of change 
of circulation term in the impulse method ts not in conformity with the behavior of the 
flow. 

The development of negative differential pressures near the central region of the 
caniber is thought to be primarilv responsible for the inception of the partial collapse 
of a parachute at high rates of deceleration. This phenomenon takes place even when 
the total drag force acting on the parachute is still positive. The sample analysis 
presented herein also shows that the negative differential pressure can cover a large 
region of the parachute and even result in negative drag. The basic idea emerging from 
the analysis reported herein is that the designs which incorporate into them the idea of 
delaving or preventing the return of the shed vortices to the canopy (€ ¢.pongsims 
management, change of deceleration historv, parachute shape, dissipation and or 
destruction of the organized wake) will be the ones which could avoid the collapse 
phenomenon. Extensive analysis and small scale experiments coupled with few 
judiciously selected field tests mav help to arrive at practically and phenomenologically 


sound parachute designs. 
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VereoNGLUSIONS 


[In this chapter the conclusions reached in the preceding sections are summarized 
and the areas in which the greatest need for further work is perceived are pointed out. 

‘A brief review has shown that two-dimensional finite difference and finite element 
methods can simulate successfully the low Revnolds number flows but serious 
impediments remain to extending them to high Reynolds number range that is of 
practical interest. The major obstacles are: the large computational requirements 
Which grow rapidly with Revnolds number, stability problems, difficulty of 
implementing boundary conditions and artificial viscosity. 

Whe <discreve volte Imocel Is scem toravold many Of the problems of the Eulerian 
finite difference methods but are subject to some of their own. Chief among these is 
the introduction of vorticity, determination of the separation points, excessive 
computation time for the convection of vortices, and the need to introduce artificial 
dissipation to bring the measured and calculated results into closer agreement. 

te meinstSSCUlOU Ol etISeStUGvedealtewiti the nilmerical simulation of a 
sinusoidally oscillating flow about a circular cylinder. 

ie mrciseicie ms Ome mO0c! "a5 §OGen combined with the bewndary laver 
calculations and the positions of the separation and stagnation points have been 
calculated as accurately as possible for a Keulegan-Carpenter number of K = 10. The 
results have accurately predicted the formation of a half Karman vortex street in the 
transverse direction. The calculated positions of the vortices were found to be in good 
agreement with those obtained experimentally. The measured and calculated drag 
force and the differential pressure distributions showed reasonably good agreement. 
The results have also shown that the effect of the backward convection of a large 
vortex over one side of the cvlinder is indeed very pronounced on all the measured and 
computed characteristics of the flow. This 1s one of the most important reasons as to 
why the Morison’s equation (see e.g., Sarpkava & Isaacson 1981) fails to represent the 
in-line force acting on the cylinder with reasonable accuracy particulariy in the range $ 
ot. = [3. 

The number of numerical parameters involved in the simulation made a complete 


Potametric study of their effects impossible. However, several tests were made in which 


ile 


the time step and convection scheme were varied. The results showed expected 
improvements as the time step was recuced and the order of convection scheme was 
increased but no undue sensitivity was observed. 

The application of the discrete vortex miodel to a particular time-dependent flow 
past a camber required an extensive study of the velocity field in the vicinity of the 
sharp edges of the camber. The results have shown that the nascent vortices can be 
introduced only at judiciously selected points. The numerical experiment predicted 
satisfactonly the evolution of the wake and the forces acting On (hewiee@ces 
Furthermore, the calculations have provided a plausible explanation for the cause of 
parachute collapse, a phenomenon which has provided the impetus for the simulation 
described herein. 

In the course of the present work, two aspects of the discrete vortex miOdel em 
identified that seem to be in greatest need for further investigation. The firstuismnme 
determination of the mobile separation points on a body without sharp edges. The 
existing methods cannot be considered satisfactory for the prediction of separation in 
unsteady laminar and or turbulent boundary lavers. The use of Pohlhausen’s method 
seemed to produce good results (as judged by the experimental data), but it 1s the most 
difficult part of the model to defend. 

The second aspect of the model! which requires further work is the discovemmuamas 
systematic and conceptually satisfactory method to reduce vorticity. It seems clear. at 
least on the basis of experiments, that there is a physical basis for dissipating vorticity. 
but 1t remains to be proved that the global error involved in doing so ts reasonable and 
in conformity with the behavior of nature. Future efforts will probably have to be 


preceded by careful measurements of the finer details of the unsteady flow field. 
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APPENDIX A 
NUMERICAL RESULTS 


- 


This section presents additional figures showing the evolution of sinusoidally 
oscillating flow about a cvlinder. It is seen that the convection of the previously shed 
forres Over ine Shoulder of the cylinder precipitates earlier separation. establishes a 
beeettcd POSition for tne mext dominant vortex, and gives rise to additional primary 


and secondary separation points. 
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Figure A.l Position of vortices, velocity and pressure distribution at T*= OES), 
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Figure A.2 Position of vortices, velocity and pressure distribution at T*= 1.100. 
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Figure A.3_ Position of vortices, velocity and pressure distribution at T*= 1.250. 
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Figure A.4 Position of vortices, velocity and pressure distribution at T* = 1.350. 
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Figure A.6 Position of vortices, velocity and pressure distribution at [*= 1.600. 
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Figure A.7 Differential pressure distributions at T*= 1.350 and 1.600. 
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Figure A.8 Position of vortices, velocity and pressure distribution at T* = 1.756. 
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Figure A.10 Position of vortices, velocity and pressure distribution at T*= 1.950. 
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Figure A.J] Position of vortices, velocity and pressure distribution at T* = 2.250. 


¥/C 


VELOCT 


PTOTAL 











me 
aes 





TANGENTIAL VELOCITY 





ROR PRESSURE 


aie 





a 


ee 


Figure A.12 Position of vortices, velocity and pressure distribution at T* = 2.300. 
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Figure A.l4 Position of vortices, velocity and pressure distribution at T* = 2.550. 
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Figure A.15 Position of vortices, velocity and pressure distribution at T*= Pat . 
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Figure A.16 Position of vortices, velocity and pressure distribution at T* = 2.925, 
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